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ABSTRACT 

An  iterative  method  for  solving  certain  boundary  value 
problems  that  occur  in  a  non-linear  theory  of  thin  elastic 
plates  and  shells  has  been  developed  [10]  and  applied  to 
obtain  numerical  solutions  of  the  von  Karman  plate  equations 
for  a  variety  of  bending  and  buckling  problems.   In  this 
report,  the  method  and  solutions  are  presented  for  circular 
plates  subjected  to  either  uniform  lateral  pressure  or  uniform 
edge  thrust  and  with  various  edge  conditions. 

The  solutions  yield  a  complete  description  of  stresses 
and  deflections  for  an  apparently  unlimited  range  of  load 
parameters.   Boundary  layer  phenomena  are  clearly  observed  and 
are  discussed.   The  Iterative  method  proves  to  be  superior  to 
the  previously  used  power  series  method  and  it  is  applicable 
to  other  non-linear  problems. 
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NON- LINEAR  BENDING  AND  BUCKLING  OF  CIRCULAR  PLATES 
Introduction. 


In  this  paper  we  consider  a  thin  circular  elastic  plate 
of  thickness   t   and  radius   R  subjected  either  to  a  uniform 
lateral  pressure,   p,   or  to  a  uniform  edge  thrust,   f.   The 
first  loading  situation  is  referred  to  as  a  bending  problem 
and  the  second  as  a  buckling  problem.   In  both  cases   only 
rotatlonally  symmetric  deformations  will  be  considered  so 
that  all  deflections  and  stresses  are  functions  of  a  single 
radial  coordinate,   r. 

It  is  well  known  that  the  classical  linear  theory  of 
plates  yields  solutions  of  these  problems  which  are  valid 
only  for  a  limited  range  of  low  loads.   For  larger  values  of 
p   and  f   it  is  necessary  to  apply  non-linear  plate  theories 
which  allow  for  finite  deflections  and  account  for  the 
interaction  of  membrane  and  bending  stresses.   Thus  we  employ 
the  von  Karman  equations  [l]  which  are  based  on  a  non-linear 
theory  of  elasticity  with  infinitesimal  strains  and  are  valid 
for  flat  plates  undergoing  small  but  finite  displacements. 
For  rotatlonally  symmetric  deformations  of  circular  plates 
these  equations  reduce  to  two  coupled  second  order  non-linear 
ordinary  differential  equations. 

The  first  extensive  investigation  of  the  bending  problem 
for  a  clamped  plate  was  given  by  Way  [2].   He  employed  power 

-  4  - 


series  to  tabulate  solutions  of  the  von  Karm^n  equations  for 
pressures  just  into  the  non- linear  range.   Subsequently, 
Friedrichs  and  Stoker  [j),^]   gave  a  complete  treatment  of  the 
buckling  problem  for  simply  supported  plates.   They  obtained 
solutions  by  means  of  perturbations,  power  series  and 
asymptotic  expansions  for  respectively,  small,  intermediate 
and  large  values  of  the  edge  thrust.   In  addition  a  boundary 
layer  effect  was  discovered  in  which  rapid  changes  of  stress 
become  concentrated  in  a  narrow  band  at  the  edge  of  the  plate 
as   f   increases.   These  methods  were  applied  by  Bromberg  [5] 
to  the  bending  problem  for  three  different  edge  conditions. 
The  boundary  layer  was  again  studied  in  detail  and  shown  to 
be  somewhat  more  complicated  than  in  the  buckling  case.   Other 
treatments,  using  the  above-mentioned  techniques,  have  been 
given  by  Chien  [6,7]  and  Bodner  [8]  for  bending  and  buckling 
problems  respectively.   In  all  of  these  studies  power  series 
solutions  could  be  obtained  only  for  limited  values  of  load. 
Green  and  Southwell  [9]  applied  relaxation  techniques  to 
finite  difference  approximations  of  a  bending  and  a  buckling 
problem  for  comparatively  low  loads.   However  the  convergence 
of  their  procedure  has  not  been  studied  and  it  is  not 
convenient  for  application  to  a  large  range  of  load  parameters. 
In  addition,  a  special  Rayleigh-Ritz  procedure  was  required 
in  treating  the  buckling  problem. 

In  this  paper  we  present  a  nxomerical  method  for  solving 
non-linear  bending  and  buckling  problems  of  circular  plates. 
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The  method  consists  of  approximating  the  von  Karman  equations 
by  corresponding  finite  difference  equations  and  then  solving 
this  non- linear  algebraic  system  by  an  "interpolated" 
iterative  procedure.   A  theoretical  study  of  the  convergence 
of  this  procedure  has  been  given  [10]  and  some  of  these  results 
are  summarized  in  the  Appendix.   This  study  and  numerous 
applications  of  the  method  indicate  that  it  yields  accurate 
solutions  for  an  unlimited  range  of  loads.   Thus  a  more 
extensive  study  of  the  development  of  boundary  layer  phenomena 
is  possible.   Furthermore,  the  method  is  applicable  in  its 
present  form  to  other  one  dimensional  non-linear  elasticity 
problems.   With  proper  modifications,  it  may  be  applicable 
to  two  dimensional  bending  and  buckling  problems. 

2.   Mathematical  Formulation. 


Let  w(r)   be  the  deflection  of  the  plate  normal  to  the 
unstrained  middle  surface,  5(r)   be  the  Airy  Stress  Function 
and  p  be  the  uniform  lateral  pressure.   Then  in  terms  of 
the  dimenslonless  variables 


(1) 


a(x)  s  -(i2(i-v2)]V2(R)  dwM  ^  p  ^[12(l-v2)]3/2  ^r^^^^ 


the  von  Karman  equations  [l],  in  polar  coordinates  for 
rotationally  symmetric  deformations,  reduce  to 
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(2a)         L  a(x)  =  -  a(x)  y{x)   -   Px^   , 


(2b)  L  y{x)    =  ^  a^(x) 


Here  we  have  introduced  the  linear  differential  operator 


^  =  ^  dx  X  dx  ^  ' 


and  in  (l)   E  is  Young's  Modulus  and  v  is  Polsson's  ratio. 

The  radial  and  circumferential  membrane  stresses,   6^ 
and  fip  respectively,  and  the  corresponding  reduced  stresses 

and 
7(x),   by 


„.^>^  ,     _,   are  given  in  terms  of  the  new  stress  function, 

■ft  — L> 


(3a)       [I,(x)   ,  Mi=l!i  (f)'  0,(r)   =  -  ^     , 


(3b)       UcM  .  iiii^  (f)'  0,M  =  -  Mxi 


■a 

The  radial  and  circumferential  bending  stresses,   cSp  and 

6p     respectively,  and  the  corresponding  reduced  stresses 

}    p  and  ^  p  are  given  in  terms  of  the  negative  slope,   a(x),   by 
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Boundary  conditions  at  the  center,   x  =  0,   and  edge 
of  the  plate,   x  =  1,   must  be  specified  to  complete  the 
formulation.   From  the  assumed  symmetry  and  regularity  at 
the  center  it  is  necessary  that 

(4)  a(0)  =  y(0)  =  0   . 

At  the  edge,  conditions  of  the  following  form  are  considered, 

(5)  a  ^^   +  bad)  =  0,   c  Mil  +  dT(l)  =  F   . 

Here  a,   b,   c  and  d  are  constants  which  define  a 

12  f ]  —  V  \       R 
particular  problem,  and  F  =  -  — ^— g — ^  (•^)   f  is  the 

dimensionless  edge  thrust.   We  consider  in  this  paper  five 

special  cases  of  (5) : 

A)  Clamped  and  radial  displacement  equal  to  zero, 
a  =  0,  b  =  1,  c  =  1,  d  =  -V,  F  =  0. 

B)  Simply  supported  and  radial  displacement  equal  to  zero, 
a  =  1,  b  =  V,  c  =  1,  d  =  -V,  F  =  0. 

C)  Simply  supported  and  radial  membrane  stress  equal  to  zero, 
a  =  1,  b  =  V,  c  =  0,  d  =  1,  F  =  0. 

D)  Simply  supported  and  radial  membrane  stress  specified, 
a  =  1,  b  =  V,  c  =  0,  d  =  1,  F  >  0. 

E)  Clsunped  and  radial  membrsine  stress  specified, 
a  =  0,  b  =  1,  c  =  0,  d  =  1,  F  >  0. 

Problems  A,   B  and  C  are  bending  problems  (?  /  O) 
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and  have  been  discussed  in  a  previous  paper  [5]-   Problems   D 
and  E  are  buckling  problems  for  which  P  =  0.   D  Is  the 
case  considered  by  Frledrichs  and   Stoker  [j>,^]   and  E  by 
Bodner  [8]. 

Equations  (2),  (4)  and  (5)  form  a  complete  statement  of 
the  problems  to  be  solved.   We  seek  solutions  for  a  wide  range 
of  the  loading  parameters  P  and  P. 

Ji.      Numerical  Formulation  and  Iterative  Solutions. 

To  obtain  niunerlcal  solutions  we  proceed  by  placing  a 
uniform  mesh  of  points 

(6)  x^  =  15x   ,   6x  =  —  ,   1  =  0,  1,  . . . ,  m   , 

on  the  dlmenslonless  radius  of  the  plate,   0  _<  x  <  1 .   The 
functions   a(x)   and  7(x)   are  then  approximated  by  "mesh 
functions"   a(x. )   and  t(x. )   which  are  defined  only  at 
points  of  the  mesh.   These  mesh  functions  are  to  be  solutions 
of  finite  difference  equivalents  of  (2),  namely: 

(Ta)  L53(x^)  =  -  a(x^)  y(x^)  -  P  x^ 

0  <  1  <  m 

(7b)  h'yi^i)   =^^2  {x^^^ 

Here  the  difference  operator,   L^.,   obtained  by  replacing 
derivatives  In  the  differential  operator  L  by  centered 
difference  quotients,  is  defined  for  any  mesh  function 
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(5x)   (^  x^  +  gSx 

X.  -  gSx  J 

At  the  endpolnts,   1=0  and  m,   the  mesh  functions  must 
satisfy  numerical  equivalents  of  the  boundary  conditions 
(3)  and  (4),  say"^ 

(8a)  a(xQ)  =  y(xq)  =  0   , 

and 

a[a(Xjjj)  -  a(x^_-j^)]  +  5x  b  a{x^)    =  0   , 

c[t(Xj^)  -  Y(Xqj_-^)]  +  5^  d  y(Xjjj)  =  6x  F   . 


(8b) 


Equations  (7)  and  (8)  are  a  system  of  2m  non-linear  algebraic 
equations  for  the  determination  of  the  2m  unknowns 
[a(x. ),  y(x.)].   The  solution  of  this  algebraic  problem  Is 
assumed  to  give  a  "good"  approximation  to  the  functions 
[a(x),  y(x)]   for  sufficiently  small   6x.   Indeed  It  can  be 
proved.  In  some  cases  [10],  that  the  mesh  functions  converge 

Conditions  (8b)  may  be  replaced  by  corresponding  ones 

centered  at  x  =  x  -  5</2. 

m     ^ 
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as   6x  — >  0  to  the  solution  of  the  continuous  problem. 

As  the  non-linear  algebraic  equations  (7)  cannot,  in 
general,  be  solved  explicitly  we  seek  iterative  solutions. 
Starting  from  some  initial  estimate,   [aQ(x.),  Yq(x. )],   of 
the  mesh  functions  for  any  value  of   P,   a  simple  iterative 
procedure  is  to  define  a  sequence  of  mesh  functions 
■(a^(Xj_),  Tj^(^i)j   ^y  ^^^   recursions 

_  _        _  2 

(9a)   Lga^(x^)  =  -  a^_i(x^)  r^.^Cx^)  -  Px. 


0  <  i  <  m 


(9b)    L^Y.(x  )  =^  a2(x,) 


Each  set  of  iterates  is  to  satisfy  (8)  at  the  endpoints,   i  =  0,  m. 
For  each  n,   equations  (9a)  now  form  a  system  of   (m-2)   linear 

equations  in  the   (m-2)   unknowns,  ^n(xj^),  i  /^  0,  m.   These 

2 
equations  may  be  solved  explicitly  in  a  simple  manner.    The 

same  is  then  true  of  the  system  (9b)  and  the  y  (x.)   may  also 

be  obtained  explicitly.   The  convergence  as  n  — >  00   of  this 

simple  iteration  procedure  has  been  exeunined  in  detail  [10] 

for  application  to  problems  A),   B)   and  C).   A  summary  of 

this  analysis  is  contained  in  the  Appendix.   For  these  cases  it 

2 

The  coefficient  matrix  of  this  system  is  of  Jacobi  or  "triple 

diagonal"  form.   Such  a  matrix  can  be  factored  into  the  product 

of  an  upper-  by  a  lower-triangular  matrix.   The  only  non-zero 

elements  of  these  factors  lie  on  the  diagonal  or  immediately 

to  its  right  or  left  respectively.   The  solution  of  the  system 

is  then  obtained  by  evaluating  two  two-term  recursions. 
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is  proved  that  the  iterations  converge  only  for  a  limited 
range  of  load  parajneters,   0  <  P  <  P',   (see  Table  I  in 
Appendix)   and  that  they  actually  diverge  for  sufficiently 
large  values  of  the  parameter  P  >  P**.   (Analogous  results 
can  be  obtained  for  the  buckling  problems   D  and  E. ) 
Since  solutions  for  loads   P  »  P**  are  of  primary  interest 
we  discard  the  simple  iterative  procedure  and  seek  one  that 
converges  for  an  unlimited  range  of  P. 

A  new  iterative  procedure  is  suggested  by  examining  a 
physical  explanation  of  the  divergence  of  the  iterations 
(9)-   Each  iterate,  a   (x.),   may  be  thought  of  as  the 
solution  for  a  plate  under  a  "non-uniform  normal  load"  in 
the  linear  theory  of  bending.   This  load  consists  of  the 
uniform  lateral  pressure  P  plus  a  non-uniform  "normal 
force"  resulting  from  the  slope  and  membrane  stress  given 

by  the  previous  iterate,   t^n-l^'^l^'  '^n- 1  ^ "'^l ^  •' '   '^^^  stress 
function,  y    (x.)^   is  a  result  of  the  deflection  determined 
by   a  (x, ) .   Now,  physically,  if  the  "non-uniform  load" 
Induced  by  a  given  iterate  is  everywhere  larger  in  absolute 
magnitude  than  the  corresponding  load  which  produced  this 
given  Iterate,  the  sequence  should  diverge.   (This  is  shown 
to  be  the  case  in  [10]  and  in  the  Appendix  when   P  >  P**.) 
Thus  it  may  be  possxole  to  eliminate  divergence  by  properly 
correcting  (9a)  for  over-estimates  in  the  next  iterate.   We 
attempt  such  a  correction  by  introducing  the  provisional 
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iterate  a*(Xj^)»   as  in  (9a),  by 


(10a)         L5a*(x^)  =  -  a^^^i^^)   Jn-i^^'i'^    -  Px^  ,   0  <  i  <  m 


and  then  Interpolating  between  it  and  the  previous  iterate  to 

obtain 

(10b)        "n^^'i^  =  ©a*(x^)  +  (l  -  Q)\.i{^i_)      • 

Here  0  is  an  interpolation  parameter  which  must  lie  in  the 
range  0  <  0  <  1.   The  iterates  y^i^^)      are  defined  again  by 
(9b)  and  both  iterates  are  to  satisfy  (8)  at  the  endpoints  of 
the  mesh.   The  solution  of  the  equations  for  the  interpolated 
iterate  case  can  be  obtained  Just  as  in  the  simple  iteration 
case  since  similar  linear  equations  need  be  solved.   It  is 
clear  that  taking  0=1   in  (lOb)  reduces  the  interpolated 
procedure  to  the  simple  iteration  (9) • 

The  convergence  of  the  interpolated  iterates  (lO),  (9b) 
and  (8),  for  proper  choices  of  0  and  initial  iterates, 
[aQ(x.),  Yq(x^)],  is  made  plausible  in  [10].  An  empirical 
method  for  determining  proper  values  of  0  is  presented  in 
the  next  section  and  a  theoretical  bound  on  0  is  given  in 
[10].  The  proposed  niimerical  procedure  is  then  well  suited 
for  application  on  automatic  computing  machines. 
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4.   Computational  Procedures. 

In  applying  the  difference  equations  (7)  It  Is  required 
that  the  mesh  spacing,   6x,   be  "sufficiently  small."  We 
determined  a  satisfactory  value  from  a  series  of  test 
calculations  with  successively  finer  meshes,  say  5x  =  0.1, 
0.05  and  0.02.   From  these  tests  it  was  observed  that  the 
.02  mesh  provided  accurate  solutions  for  the  ranges  of  loads 
considered. 

Application  of  the  iteration  scheme  (lO) ,  i9h)    to  solve 
the  difference  equations  (7)  requires  procedures  for  determining 
an  initial  iterate,  a  "best"  value  for  ©  and  a  criterion  for 
convergence  of  the  iterates.   For  the  latter  we  use 

(11)   max  (|a^(Xi)-an-l^^i)l  '  '^n^^i^" Vl^'^i^  I  J  <  ^°'^  ' 

and  iterate  until  this  condition  is  satisfied.   It  is  found 
that  the  number  of  iterations  required,  starting  from  a  given 
initial  Iterate,  depends  strongly  upon  the  value  of  0  used 
in  (lOb).   A  typical  graph  of  this  dependence  of  the  number 
of  iterations  for  convergence,   N(©),   on  Q     for  a  given 
value  of  the  load  parameter  is  shown  in  Fig.  1.   The  "best" 
value  of  ©,   is  that  value,   ©q^.^   for  which  N(©)   is  a 
minimum.   From  such  graphs  for  a  sequence  of  values  of  the 
load  parameter,  say  P,   a  corresponding  sequence  of  values 
of   ©on^^^   "^^  obtained.   From  a  smooth  curve  through  the 
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points   (P,0   (P))   good  estimates  of  ©    for  an  entire 
range  of  load  may  be  determined  (see  Fig.  2).   It  was  found 
that  0    does  not  change  appreciably  for  significant 
changes  In  6x  and  thus  curves  of   ©on^^^   were  determined 
with  few  mesh  points  and  little  computing  time.   A  graph  of 
0   (?)   for  problem  A)  Is  shown  In  Fig.  2,    together  with  a 
theoretical  upper  bound  determined  In  [10].   Similar  graphs 
have  been  obtained  for  the  other  problems. 

The  Initial  Iterate  for  a  given  load  P  In  any  of  the 
bending  problems,  A),  B)  or  C) ,  is  taken  to  be  the  final 
("converged")  iterate  for  a  slightly  smaller  value  of  load, 
that  is 

aQ(P,x^)  =  "^(P  -  5P,x^)   . 

In  this  manner  solutions  were  obtained,  with  the  aid  of  the 
previously  estimated  ©oqCP)*  ^o^  ^^  apparently  unlimited 
range  of  load. 

The  buckling  problems,  D)  and  E) ,  have  non-unique 
solutions  and  so  the  above  procedure  must  be  applied  with 
care.   It  has  been  shown  by  Frledrlchs  and  Stoker  [3>^] 
that  for  F  <  Fq   (where  Fq  is  the  first  eigenvalue  of  the 
linear  buckling  problem)  the  only  solution  is  the  trivial  one: 
a(x)  =  0  ,  y(x)  =  Fx.   For  F  >  F^  two  solutions  [that  differ 
from  each  other  only  in  the  sign  of  a(x) ]   in  addition  to 
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the  trivial  one  are  possible.   If  the  trivial  solution  is 
uaed   aa  an  initial  iterate  it  might  then  result  as  the 
solution  for  all  F  and  the  non-trivial  solutions  would 
not  De  detected.   To  avoid  this  possibility  we  use  the  non- 
zero initial  iterate  ciqCx,)  =  x,   for  some  F  <  Fq  and  the 
iterations  then  converge  to  the  trivial  solution.   This 
procedure  is  continued  for  a  sequence  of  increasing  values 
of   F  until  a  non-trivial  solution  is  obtained.   This  non- 
trivial  solution  is  then  used  as  the  initial  iterate  for  the 
next  value  of   F   and  the  procedure  is  continued  as  in  the 
bending  problems.   In  this  manner  a  knowledge  of  the  value  oi' 
F^      is  not  required  and  it  can  in  fact  be  estimated  from  the 
computations.   It  is  also  found  that  the  value  of  ©^  (F) 
for  obtaining  trivial  solutions  when  F  <  F^  is  greater  than 
one  and  extrapolation  is  automatically  used  in  (ICb) .   Near 

F  =  F^,   0   (F)  »  1   and  as  the  non-trivial  solutions  are 
0    op^  ' 

oDtained  for  increasing  F,   0  _(F)   decreases.   Thus  it 
would  seem  that  the  curve  of   ®or|(F)   ^^^  ^  bifurcation,  as 
do  the  solutions,  at   F„. 

When  the  iterations  converge  for  a  given  load  the 
stresses  and  deflections  can  be  computed.   For  this  purpose 
we  employ  obvious  finite  difference  equivalents  of  equations 
(l)  and  CJ)  and  evaluate  them  using  the  final  iterates, 
['aj,(x.),  ^^(x.  )  ] .   Such  results  are  presented  and  discussed 
in  the  next  section. 
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3.   Discussion  of  Results.  -^ 

Representative  graphs  of  the  reduced  stresses  and 
deflections  as  functions  of  x  for  Increasing  sequences  of 
load  are  presented  In  Fig's,  (^a-e)  for  bending  problem  C), 
and  In  Fig's.  (4a-e)  for  buckling  problem  E) .   Similar  results 
have  been  obtained  for  the  other  three  problems  and  those  of 
problem  A)  are  presented  In  [10]. 

In  Fig.  3b,  the  reduced  slope,   a,   for  problem  C)  Is 
shown  to  have  a  maximum  near  the  edge.   Although  It  Is  not 
obvious  from  the  curves,  the  data  Indicates  that  the  location 
of  this  maximum  moves  towards  the  edge  as   P  Increases.   The 
reduced  circumferential  membrane  stress  (Fig.  5c)  remains 
tensile  and  rather  flat  In  the  center  of  the  plate  and  changes 
rapidly  to  compressive  with  Its  maximum  at  the  edge.   The  point 
of  zero  stress  moves  toward  the  edge  with  Increasing  P.   This 
sharp  transition  of  the  stress  near  the  edge  from  tensile  to 
compressive  Is  characteristic  of  the  "boundary  layer"  behavior 
first  discovered  by  Friedrichs  and  Stoker  [5,4].   The  radial 
membrane  stress  (not  shown)  has  a  graph  similar  to  the  lateral 
deflection  W  of  Fig.  3a  with  a  maximum  at  the  center.   It 

All  the  numerical  calculations  were  performed  on  the  Atomic 
Energy  Commission  UNIVAC  at  NYU.   The  authors  wish  to  express 
their  appreciation  to  M.  Weissner,  I.  Krongold  and  S.  Weinberg 
for  their  aid  in  coding  these  problems  for  the  machine. 

4 

In  all  calculations  Poisson's  ratio,   v  =  .3. 
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can  be  observed  from  these  graphs  that 

(12)         maxl^^l  >  maxlTZj^l   . 

For  low  values  of  P,   the  bending  stresses,  Fig's. 
(3d,e),  have  their  maxima  at  the  center  of  the  plate.   As 
P   increases  these  points  of  maximum  stress  move  toward  the 
edge  thus  again  exhibiting  a  boundary  layer  behavior.   In 
addition  we  observe  that 

(15)  max|y~g|    >  max|2_pj|    >  max|/_|-,|       , 

where  the  equality  holds  only  when  the  maximum  bending  stresses 

are  at  the  center.   Hence  the  dominant  stress  is  the  circumferential 

membrane  stress  at  the  edge. 

The  curves  for  bending  problem  A)  are  presented  in  [10]. 
Here  the  graphs  of   a  clearly  show  that  the  location  of  the 
maxima  move  towards  the  edge  as   P   increases.   Both  membrane 
stresses  remain  tensile  with  their  maxima  at  the  center  and 
hence  their  maxima  are  equal.   The  bending  stresses  are  tensile 
and  almost  constant  in  the  central  portion  of  the  plate  and 
rapidly  change  to  compressive  and  maxlmxim  at  the  edge.   Again 
the  points  of  zero  stress  approach  the  edge  with  increasing 
P.   The  dominant  stress  in  this  problem  is  the  radial  bending 
stress  at  the  edge  of  the  plate. 

For  problem  B)  the  graphs  of  the  slope  and  bending 
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stresses  are  similar  to  those  of  problem  C)  and  the  membrane 
stresses  are  similar  to  those  of  problem  A) .   The  maxima  of 
all  the  stresses  are  related  as  In  (13)- 

The  reduced  lateral  deflection  and  slope  for  buckling 
problem  E)  (see  Fig's.  4a, b)  flatten  in  the  central  portion 
of  the  plate  as  F  increases.   Thus  in  this  central  region 
the  slope  approaches  zero  while  in  the  bending  problems  the 
slope  increases  as  the  load  parameter  increases.   For  F 
slightly  larger  than  the  critical  value  the  membrane  stresses 
(see  Fig.  4c)  are  tensile  in  the  central  region  and  flatten 
in  this  region  as   P   increases.   These  stresses  rapidly 
decline  and  become  compressive  reaching  a  maximum  at  the 
edge  (for  the  radial  stress)  or  near  the  edge  (for  the 
circumferential  stress).   The  maxima  of  the  membrane  stresses 
are  related  as  in  (12) .   The  bending  stresses  are  tensile  in 
the  centraj.  region  reaching  a  tensile  maximum  near  the  edge 
and  then  rapidly  decreasing  to  compressive  with  the  maximum 
at  the  edge,  see  Fig's.  (4d,e).   These  stresses  approach  zero 
with  increasing  F,   in  the  central  region,  and  both  the 
locations  of  tensile  maximum  and  zero  stress  move  towards  the 
edge.   Although 


max  \YZr\    =  max  |XI?' 


-R' 
equations  (3)  show  that 
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max  1 6^:j|  >  max  j  6^1 

Hence  max  |6:j|   Is  the  dominant  stress. 

The  graphs  of  lateral  deflection   (W) ,   slope   (a), 
radial  membrane  stress  (}    „)   and  circumferential  membrane 
stress  {)     „)   for  problem  D)  are  very  similar  to  those  for 
problem  E) .   However  the  bending  stresses  are  always  tensile 
with  the  maximum  moving  toward  the  edge  as  F  increases. 
The  dominant  stress  is  the  circumferential  membrane  stress 
at  the  edge.   These  results  for  problem  D)  are  In  good 
agreement  with  those  of  Friedrichs  and  Stoker  [^]. 

Figures  5a- e  are  graphs  of  the  load  parameter  vs.  the 
lateral  deflection  at  the  center  and  the  maximum  membrane  and 
bending  stresses  for  each  of  the  five  problems. 

To  compare  methods  power  series  solutions  [2,5]  were 
computed  on  the  UNIVAC  for  problem  A) .   All  of  these 
solutions  are  in  excellent  agreement  with  those  obtained 
by  the  iteration  method  presented  in  this  paper.   However 
the  power  series  method  has  two  distinct  disadvantages : 
the  computing  time  it  requires  far  exceeds  that  for  the 
present  method;  power  series  solutions  could  be  obtained 
only  for  P  <  2000.   The  present  method  easily  yields 
solutions  for  loads  as  high  as   P  =  7000  and  it  seems 
unlimited  in  this  respect. 

In  computing  the  ultimate  load  of  a  buckled  rectangular 
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plate  it  is  customary  [11]  to  assxime  that  In  the  post-buckling 
region  the  entire  load  is  resisted  by  two  narrow  strips 
adjacent  to  the  edges  parallel  to  the  loading.   In  these 
strips  it  is  further  assumed  that  the  stress  normal  to  the 
width,  i.e.  in  the  direction  of  loading,  is  constant  and 
equal  to  the  edge  stress.   The  width  of  these  strips,  which 
must  be  determined,  is  commonly  referred  to  as  the  "effective 
width, "  The  concept  of  effective  width  is  clearly  related 
to  the  boundary  layer  and  its  characteristics  in  the  neigh- 
borhood of  the  edge.   In  fact  some  measure  of  the  "boundary 
layer  thickness"  may  be  considered  to  be  the  effective 
width  and  conversely. 

From  the  results  for  the  circular  plate  it  is  possible 
to  gain  some  Insight  into  appropriate  measures  of  the 
boundary  layer  thickness.   From  analogy  with  the  rectangular 
plate  we  may  assume  in  this  case  that  the  stresses  perpendi- 
cular to  the  boundary  layer,  the  circumferential  stresses, 
are  of  primary  importance.   For  example  in  problem  E) , 
Fig's,  (4c  and  e),  the  distance  from  the  edge  to  the  point 
of  zero  circumferential  stress,  which  decreases  with 
increasing  F,   may  be  taken  as  a  measure  of  the  boundary 
layer  thickness.   This  distance  is  plotted  as  a  function  of 
F  in  Fig.  6.   From  these  graphs  empirical  formulae  for  the 
variation  of  the  thickness  with  F  can  be  obtained. 
Alternate  definitions  determined  by  other  characteristics  of 
the  boundary  layer  could  have  been  employed. 
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APPENDIX:  Theoretical  Convergence  Properties. 

In  this  section  a  summary  of  the  convergence  properties 
of  the  simple  Iteration  procedure  Is  presented.   Only  the 
bending  problem  A  will  be  considered.   All  results  obtained 
for  this  case  have  exact  analogies  for  the  other  bending 
problems,   B  and  C,   and  a  similar  analysis  can  be  carried 
out  for  the  buckling  problems   D  and  E. 

The  formulation  given  by   (2),   (4)   and   (5)A)   can  be 
replaced  by  the  equivalent  Integral  equations 


1 
a(x)  =  f(x)  -/  g^{x,i)a(Oyi0^i     . 


(A.l) 


where 


o 

1 

I 


7(x)  =  1/  g2(x,0a^(0d^  i 


(A. 2)    f(x)  =|x(l-x^)  =/  gi(x,OP|^de   ; 

o 


and 


(A.5) 


g2(x,0  =  -  2  ^  '  H  =  1  +  V 


-  24 


From  the  differential  equations  and  the  fact  that  gj^(x,0 
and  goCxjO  <  0  ^or  0  <  x,^  <  1  we  may  conclude  that 

0  <  a(x)  <  f (x) 

(A. 4)  ^ 

0  >  7(x)  >  g(x)  =  ^Jg^{x,i)r^ii)d   =  -i^f^ 

[{k   +  i)-x2(6-4x2+x'^)] 

These  are  bounds  on  the  exact  solution  and  are  valid  for  all 
loads. 

The  n-th  Iterate  of  the  numerical  procedure  (9)  is  a 
numerical  approximation  of  the  n-th  iterate  of  the  analytical 
iteration  scheme  defined  by:   a^Cx)  =  0  , 


(A.5) 


L  a^(x)  =  -^n-l^^^'^n-l^^^  "^^ 


L  y^M    =  2  °^n^^^ 


where  each  iterate  is  to  satisfy  the  boundary  conditions  (5) A). 
Since  the  solutions  of  these  equations  are  sufficiently  smooth 
it  can  be  proven  that  the  numerical  solutions  of  (9)  converge 
as   6x  — >  0  to  the  solution  of  (A.5)-   Thus  the  convergence 
of  the  numerical  iterations  (as  n  — >  oo  )  is  equivalent  to 
the  convergence  of  the  analytic  iterations. 

The  solution  of  (A.5)  is,  in  terms  of  the  Green's 
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functions  (A. 3), 

1 

o 
(A. 6) 

o 


If   a2(x)  >  0  for  0  ^  X  <  1  we  may  then  prove  by  induction 
that 


0  =  Oq  <  CX2  < <  °^^  < <   °^2n+l  ^ ^  a^(x)=f  (x) 

(A.7)  fO<x<l 

0  =  Yq  >  72    >•  •  •  •>  Ygn  > >  72n+i    > >  'Yi(x)=g(x)^ 

The  bounds  of  these  sequences  are  Just  the  bounds  (A. 4)  on 
the  exact  solution.   E^r  a  simple  calculation  It  can  be  shown 
that  a2(x)  >  0  If  P  <  P*  ,   where  P*  can  be  explicitly 
evaluated  (see  Table  I) . 

If  a2(x)  <  0  for  all   0  <  x  <  1   the  Iterations  (A. 6) 
must  diverge,  forming  the  sequences 


•  •  •  •  <  °^2n  ^ <  ^2  <  Oq  <  a^  < <  a^^^^   <• 

(A.8)  yo^  'Yl  ^ ^T'2n+1  ^• 


The  previous  calculation  also  shows  that  ap(x)  <  0  if 


-  26  - 


P  >  p**   (see  Table  I  for  p**) . 

It  is  possible  that  for  P  <  P*   the  Iterates,   a^(x), 
will  not  converge,  I.e.  the  even  and  odd  subsequences  In 
(A. 7)  could  converge  to  different  limiting  functions. 
However,  we  can  prove  that  there  Is  some  P'  <  P*  such  that 
they  converge  for  all   P  <  P'.   For  this  proof  consider 
I  a  -,(x)  -  a  (x)|   and  use  the  bounds  f(x)   and  g(x)   of 
(A. 7)  to  derive 


(A. 9)         max|a^^-^(x)  -  a^(x)|  <  k"  max|f(x)|   , 

X  X 

where 

P  2 
(A. 10)         K  =  irr)      max  p,(x)   , 

°   0<x<l   -^ 

and  p^(x)   Is  a  non-negative  polynomial,  independent  of 
P,   in  0  <  X  <  1.   Taking 

_1 

P^(x)]' 
0<x<l 


(A. 11)        P'  =  8[  max  p^(x)]     <  P*   , 


we  have  K  <  1  for  all   P  <  P'   and  hence  by   (A. 9)  the 
sequences  (A. 7)  converge  as  n  — >  oo  . 

The  values  of  P',   P*  and  P**  are  given  in 
Table  I  for  problems   A,   B  and  C   A  more  detailed  estimate 
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Prob.  A 

Prob.  B 

Prob.  C 

P' 

68.5 

9.9 

24.0 

p* 

1250 

14.9 

28.3 

p»* 

197.0 

15.6 
Table  I 

42.6 

of  K  in  deriving  (A. 9)  would  raise  the  values  of  P' 
above  those  presented  in  the  table. 
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Typical   Graph  of   Number  of   Iterations    for    Convergence. 
Problem   A,     P=  300,     Sx  =  .05,     a  Ax.)  -   const. 
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Fig.  3a'     Reduced  Lateral  Deflection,   Problem  C 
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Fig.  3b"-    Reduced  Slope,     Problem  C. 


Fig.  3c'    Reduced   Circumferential  Membrane    Stress, 
Problem  C. 


Fig.   3d:    Reduced    Radial  Bending    Stress, 
Problem   C. 


Fig.  3e-    Reduced  Circumferential  Bending  Stress, 
Problem  C. 
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Fig.  4a-    Reduced  Lateral  Deflection,    Problem  E 
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Fig.  4b:    Reduced    Slope,    Problem    E. 


Fig.  4c-     Reduced  Circumferential  Membrane   Stress, 
Problem    E. 


Fig.4d:      Reduced    Radiol  Bending    Stress, 
Problem    E  , 
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Stress,    Problem   E. 
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Fig.  5b-   Maximum    Stresses    and  Deflection    at 
the  Center.    Problem  B. 
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Fig.  6:    Variation   of  Boundary  Layer    Thickness,  Problem   E 
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